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Abstract 
T3 ■ 

We study numerically the crossover between organized and disor- 
ganized states of three non-equilibrium systems: the Poisson/coalesce 
random walk (PCRW), a one-dimensional spin system and a quasi 
one-dimensional lattice gas. In all cases, we describe this crossover 
in terms of the average spacing between particles/domain borders 
{S(t)) and the spacing distribution functions p^(s). The nature of 
. the crossover is not the same for all systems, however, we found that 

^sf \ for all systems the nearest neighbor distribution p'°'(s) is well fitted 

^sO • by the Berry-Robnik model. The destruction of the level repulsion in 

the crossover between organized an disorganized states is present in 
| all systems. Additionally, we found that the correlations between do- 

mains in the gas and spin systems are not strong and can be neglected 
in a first approximation but for the PCRW the correlations between 
rjj '■ particles must be taken into account. To find p( n \s) with n > 1, we 

propose two different analytical models based on the Berry-Robnik 
model. Our models give us a good approximation for the statistical 
behavior of these systems in their crossover and allow us to quantify 
the degree of order /disorder of the system. 
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1 Introduction 

Many one-dimensional non-equilibrium systems exhibit a crossover between 
disorganized and organized states. This crossover usually depends on the 
value of one or more parameters which set the system in one of these states. 
Our main objective is to study the statistical behavior of three non-equilibrium 
systems in their crossovers between organized and disorganized states. The 
first system is the Poisson/coalesce random walk (PCRW) where the parti- 
cles describe independent random walks and when two particles meet they 
coalesce with probability k, otherwise, they interchange their positions. The 
second system is a quasi one- dimensional gas, where the particles interact 
only by volume exclusion in presence of an external field. The last system 
is a one- dimensional spin lattice, where the particles interact by a coupling 
force J in presence of an external driving field. 

We propose two analytical models for the spacing distribution functions 
of these systems and compare them with the numerical results from the 
simulation. Our analytical models are based on the Berry-Robnik model 
introduced in Ref . pQ . This model is used in quantum systems which are nei- 
ther integrable nor fully chaotic in order to find an analytical approximation 
to the nearest neighbor distribution of the energy levels of the system, see 
Refs. [IJ[2j E]. The Berry-Robnik model depends on one parameter which con- 
trols the crossover. The crossover is described through the nearest neighbor 
distribution p^'(s), which gives us the probability that the distance between 
two consecutive levels is s. We complement the model proposed in Ref. [1] 
calculating analytically the higher order spacing distribution functions p( n \s) 
for n > and the pair correlation function g(r) by two methods. 

This paper is organized as follows. Because of its importance, in the 
second section we explain in detail the Berry-Robnik model. In third, fourth 
and fifth sections, we test our analytical models with three non-systems: 
the Poisson/coalesce random walk, the quasi one-dimensional gas and one- 
dimensional spin lattice respectively. These systems are also explained in 
those sections. 
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2 Berry-Robnik model as a non-equilibrium 

system 

We consider a continuous one-dimensional ring with two kinds of particles A 
and B and normalized densities pA and ps respectively. The A particles are 
subject to the reaction A + A — > A, i.e., they are coalescing random walkers. 
Particles B describe independent random walks and they are subject to the 
reaction B — > 0. The amount of B particles that disappear is calculated in 
such way that the number of A particles, Na, and the number of B particles, 
Nb, satisfy the relation Nb/Na = constant at each time step. This is the 
only interaction between both species of particles. 

For this system, P^ n '(S,t) is the probability that the distance between 
two particles is S at time t under the condition that between these particles 
there are n additional particles. When the average distance between neighbor 
particles is much smaller than the size of the lattice, the system exhibits a 
dynamical scaling for all values of pa and ps- In this regime the spacing 
distributions can be scaled by using the variable change s = S(t)/ (S(t)). 
Then, we obtain the time independent spacing distributions pW(s). 

It is known that in the scaling regime Na = L/ (27iDt) 1 ^ 2 where D is 
the diffusion constant and L is the size of the lattice, see Ref. [I] for more 
information about the coalescing random walk (CRW). It is straightforward 
to show that Nb = NaPb/{^ — Pb), then, the average length between nearest 
neighbors satisfy (S(t)} = (1 - p B )(2rrDt) 1 / 2 . 

The interaction between species warrants that the quotient between nor- 
malized densities pA and p B is constant for all time, i.e., the rate of disap- 
pearance of particles B is proportional to the one for A particles and the 
proportionality constant is Pb/pa- Taking this into account, in the scal- 
ing regime, the system is the uncorrelated superposition of two independent 
systems of particles A and B with constant normalized densities. Our ob- 
jective is to calculate the nearest neighbor distribution p^°\s) of the whole 
system. This can be done easily because we have the superposition of two 
uncorrelated distribution functions pf^ (s) which satisfies 




(1) 



and 




(2) 
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where % = AorB and pi is the density of particles corresponding to the 
systems A or B. As is natural, the normalized densities satisfy p A + p B = 1- 
Let E^(s) be the probability to choose randomly an empty segment of 
length s. This probability is related to p(°\s). In order to prove that, we 
consider the probability that there are no particles in the interval q to q + r, 
given that there is a particle at q. This probability is given by J r °° p^ (x) dx. 

In our case we have the uncorrelated superposition of two kinds of parti- 
cles, then we can write 



p ( >(x)dx = p A q B (r) p\ '(y)dy + p B q A (r) P B (y)dy, (3) 

J r J r 

where, for example, p A is the probability that the particle in position q 
belongs to the A system, {x)dx and g_e(r) are the probabilities that 

there are no particles of the systems A and B in the interval, respectively. 
In order to find we choose a particle randomly, the probability that 

the chosen point q lies within a gap of length a to a + da is proportional 
to apg\a), to normalize this equation we use the fact that (s B ) = l/p B , 
obtaining the probability distribution p B daap^\a). Now, the probability 
that the distance to the next particle is r, given that the point is in the gap 
of length a, is zero if a < r and 1/cr if 1 < r < a. The probability of not 
having a particle of the system B in an interval of length r is (1 — r/a)6(a — r) 
{0 is the unit step function). 

Thus, the unconditional probability that the distance until the next par- 
ticle is r is given by 

/oo 
p { B \(y)((r -r)da. (4) 

The same argument can be done for the second term of equation ([3]). Thus 
Eq. (E]) takes the form 

) /»oo /*oo 

p (S) \x)dx = pa P { A{y)dyp B I pf' \<t){<t - r)da 

J r J r 

POO POO 

+ Pb P B \y)dyp A I p { A \a)(a -r)do. (5) 

J r J r 

The probability E^(s) is the probability that r > 5, then, integrating Eq. (El) 
over r we find 

/oo /*oo 
dr ^ dxp {Q \x)=p A E { l\s)p B E^\s) ) (6) 
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where we define 



DC 



E?\s)= dapf\a)(a - s) = da dxpf\x). (7) 

J s J s J a 

From equations (jSJ) and ([7j), the spacing distribution function p^(s), 
which results of the mix of systems A and B is given by 



d 2 E^[ 



s 



p( o, (s) = (8) 

The fact that E^(s) is given by the product of E t (s) functions is nat- 
ural, because we used two statistical independent sequences. It is easy to 
prove that the E^^s) functions satisfy 

^ (1) (0) = ^, (9) 
Pi 

and 

dE^(0) , . 

— = -i- (10) 
ds 

From the above equations it follows that p^ (0) is given by 

P {0 \0) = l + PA ( P f(0) - Pa) + Pb (p^CO) - PB ) . (11) 

Then, even if p^(0) = and Pg\o) = it can happen that p(°)(0) ^ 0. In 
our particular case, we superpose one Wigner distribution 

P {0) (s,p A ) = ^p 2 A seM-lpW] (12) 
with density pA and one Poisson distribution 

p (0) (s,p B ) = PBexp[-p B s] (13) 
with density p B , then, we have 

pV>(s,q) = e-* s (Verfc (^(1 - q)s^j + (2q(l - q) + |(1 - q) 3 s) e -f(W)V ) |(14 ) 

where p B = q is the density of Poisson sequence, pA = 1 — <? is the one 
for the Wigner sequence and erfc(z) is the complementary Gaussian error 
function. Note that Eq. (fUjl reduces to the Poisson distribution for q = 1, 
and for q = it reduces to the Wigner distribution. This result is well-know 
in random matrices theory and corresponds to the Berry-Robnik model for 
crossover between chaotic and non-chaotic behavior in quantum systems, see 
Ref. p. 
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2.1 Higher spacing distribution functions 

In Refs. [3 16] the authors show that the nearest neighbor distribution is not 
enough to describe the complete statistical behavior of a non-equilibrium sys- 
tem, because several different systems could share the same nearest neigh- 
bor distribution. Because of that, we generalize the Berry- Rob nik model 
for higher spacing distribution functions. Let E^ n '(xx,yx, ■ ■ ■ ,x n ,y n ) be the 
joint probability that the intervals [xi, yi] (i = 1,2, ••• , n) are empty. The 
intervals are non overlapping and ordered xi < y\ < • • • < x n < y n . Be- 
cause of the independent superposition nature of the Berry-Robnik model, 
E^ n '(xi,yx, ■ ■ ■ , x n , y n ) can be written as 

£ (n) (xi,?/i, • ■ ■ ,x n ,y n ) = p'%E ( x ) (x 1 ,y 1 , ■ ■ ■ ,x n ,y n )p n B E i z\x 1 ,y 1 , ■ ■ ■ ,x n ,y n ), 

(15) 

where E^\x 1: y 1: ■ ■ ■ ,x n , y n ) and E%\xi, y 1: ■ ■ ■ ,x n , y n ) are the joint prob- 
abilities for particles A and B respectively. The authors of Ref. [TJ, showed 
that spacing distribution functions p( n \s) can be calculated from Eq. (fi3|) . 
let us summarize their most important result. Let p^ n \s) be the probability 
that given one particle its (n + l)-th neighbor is at a distance s. From its 
definition p( n \s) is given by 

P in) (s)= [ a n+2 \0, yi ,--- ,y n ,s)d yi ...dy n , (16) 

J 0<j/i<-<j/„<s 



with 

J n) (x ir -- ,x n ) 



(„)/_ , d nE{n X) {xx,y u --- ,a;„-i,j/„_i) 



dxi ■ ■■dx n - 1 dy n ^ 1 



• (17) 

yi=x2,--- ,y n -i=x n 



On the other hand, it is well-know that for CRW, E^ (xi, Vi, • • • , x n , y n ) 
is given by [I] 

E^\xi, yi, ■ ■ ■ ,x n ,y n ) = a p E^\zi tP , Z2, p ) • ■ • E^\z2 n -i,p, z 2n,p), (18) 



where zi tP , Z2, p , Z2 n , P symbolize an ordered permutation, p, of the variables 
x u y x , ...,x n ,y n , such that 

z l,p < z 2,pi z 3,p < z 4,p, ' ' ' i z 2n-l,p < z 2n,p, (19) 
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z l,p < z 3,p < z 5,p < ■ ■ ■ < Z2n-l,p- 

(20) 

In Eq. ffTgj) cr p is the signature of the permutation, i.e., a p =l for even per- 
mutations and Op = — 1 for odd permutations. 

The function yi) is the probability that from x\ to yi the lattice is 

empty. Then it is possible generate the complete solution for the CRW from 
E^\xi,yi), which is given by the solution of the diffusion equation under 
the suitable boundary conditions (see Ref. [1]). In fact, the exact expression 
for this function is 

E^izuy!) = -^erfc - Xl )j , (21) 

for additional information see Refs. [I] and [8]. The particles B describe 
independent random walks, then, we have 

E ( i\x 1 ,y 1 ) = —e-^-^\ (22) 
Pb 

and 

n 

Ef{x u 2/i, • ■ ■ , x n , y n ) = J] _ e -PB(w-«i). (23) 

i=i Ps 

From equations (TT5l). (fl8l), (12T1) and (1231) we can calculate E^ n ^{x\^ y±, ■ ■ ■ , x n , y n 
for the whole system and then, we can use (fl6|) and (ITTj) to calculate p( n '(s). 
By using this formalism it is also possible to calculate the n-point correlation 
function 

p<»W- )Xb) = ( _i r ______ ,*»,y»)|^,..., fc= ^, 

(24) 

in fact, the pair correlation function is given by 



g (r) = l-(l- qfe-^ 1 "^ 2 + ^e^ 1 ^ {1 - g) 3 rerfc ^(1 - q)r j . 

(25) 

The above equation takes the form g(r) = 1 for q = 1. For q = we recover 
the pair correlation function of the CRW, see Ref [_]. 



7 



3 A simple model of crossover between or- 
ganized and disorganized states: The Pois- 
son/coalesce random walk 

Consider a one- dimensional ring with L sites and n p particles, then, the parti- 
cle density is given by p = n p /L. The particles describe independent random 
walks and when two particles meet they coalesce with probability k, other- 
wise, they interchange their positions. This system was studied previously 
in Refs. [HI HD]- The algorithm used in the simulation is the following: 

1. n p particles are randomly inserted in a L sites lattice. 

2. A particle is chosen at random. 

3. The particle can move to the left or to the right with same probability. 
If the particle is moved to an occupied site, the particles coalesce with 
probability k or they interchange their positions with probability 1 — k. 

4. In a time unit all particles are moved. 

In the limit t ^> 1, finite systems reach a non-equilibrium steady state 
(NESS), where there is only one particle which executes a simple random 
walk and (S(t)) = L. As we will see soon, for an infinite size system in the 
same limit, the average length between nearest neighbor particles grows as 
t 1 / 2 for k > and the system is statistically equivalent to the CRW. 

It is possible to derive an approximate analytical solution by using the 
inter-particle distribution function method (IPDF), see Refs. [HE]. Let E n 
be the probability to find an empty segment of length n in the lattice. The 
master equation for E n can be written as 

= ^ (E n+1 (t) - 2E n (t) + E n _ x {t)) + — (k - 1) P(^?..), 

(26) 

with boundary conditions E (t) = 1 and E^t) = 0. The probability 

n-1 

P(o~^"o««) cannot be written in terms of E n {t). However in Refs. pQ, [9] 
the authors propose an approximation for this probability 



PM ~ (2?) 
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which can be written in terms of E n (t) as 

. (l-2E 1 +E 2 )(E n _ l -E n 



I -Ex 



(28) 



This approximation allows us to calculate the concentration of particles 
c(t) = 1/ (S(t)). From Eqs. (1261) and (1281) and summing over the index n 
it is easy to find 



^-^)-^))-^-*) (1 -f-y M) ^)- 

n=l v ' 

(29) 

Taking into account that d t E\(t) = 2Dk(l — 2Ei + E^/Ax 2 and making the 
approximation -EnCO ~ (2/ tt) / ( 1 — Ei(t)), in Ref [9] it is found that 

2d f 1 \_ c(T) + (l-k)dc(r 1 



71 dr \c(r) J kc(r) dr 

where r = 2Dt/Ax 2 and c(r) = 1 — -B'i(r). Equation fl30l) can be integrated 
obtaining 

c(r) = 2C ° fc (31) 

c 2 (A; - 1)tt + x/c^-Sfc + c (fc - l)vr) 2 + 4c^ 2 7rr 

with Co = c(0). Note that for r — > 00, we have c(r) oc r^ 1 / 2 as we mentioned 
above. In figure [H we show the behavior of (S(t)) = l/c(r) for different 
values of k, this result was first shown in Ref. [9]. The agreement between 
Eq. (|3"T1) and the simulation is very good. Because of the finite size effects 
in the simulation, for L = 500, the system does not reach the asymptotic 
exponent (3 = 1/2 ((S(t)) oc t") but for L = 10000 we can see this regime for 
k = 1, k = 0.25, k = 0.02 and k = 0.01. For low values of k, the time that 
the system remains disorganized with (S(t)) almost constant, is larger than 
for large values of k. 

The nearest neighbor distribution (s) evolves in the following way. The 
system starts in an disorganized state in such way that p^(s) is described 
by the Poisson distribution, then, for k the system evolves and p(°\s) 
is deformed continuously until it reaches the Wigner distribution, for large 
systems. For small systems, the finite size effects appear before this regime is 
attained, i.e., small systems reaches the NESS without reaching the scaling 
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Figure 1: Behavior of (S) as function of t. In the simulation we took two 
lattices with L = 10000 (filled symbols) and L = 500 (empty symbols) over 
500 realizations. In both cases the initial density was 1/4. The continuous 
line represents Eq. (1HT1) . 
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regime. In the k = case, the system remains disorganized for all values of 
t. 

In the continuum limit, Eq. fl26|) takes the form 



dE^(x,t) 
dt 



2D 



d 2 E^(x,t) 



+ 2D(k- 1) 



d 2 EW{x,t) 

dE<V(x,t) dx~ 2 



x=0 



(32) 



dx 2 



Ox dEW(x,t) 
dx 



x=0 



with boundary conditions E^(0,t) = 1 and E^>(oo,t) = 0. This equation 
can be solved easily for k = 1. For arbitrary values of k, Eq. (|32|) cannot be 
solved exactly, however some approximate expressions were found in Refs. [U 
[9j [10], but they involve self-consistent forms which are difficult to handle. 
Motivated by this fact, we propose to use the Berry- Robnik model to find an 
approximate analytical solution to the statistical behavior of this system. 

Equation (|14p provides a good fit for the crossover in this reaction- 
diffusion model during its time evolution as we can see in figure |2j For 
the nearest neighbor distribution the fit is almost perfect for low and high 
values of s, but for intermediate values we can see little differences in the 
interval 0.25 < s < 1 for intermediate values of k. In this figure all data 
was taken at different times with k = 0.05 over 20000 realizations, the initial 
density was 0.1. The fit parameter is q. We found the appropriate value of 
q by using the numerical results for p(°\s), equation ()14p and the minimum 
square criteria 0. 

The higher spacing distributions have been calculated with two methods. 
In the first one, we use our extension for the original Berry-Robnik Eq. (1151) 
with equation (TT6T) . From now on, we call it the generalized Berry-Robnik 
(GBR) model. For the second one, we use equation (fT4|) with the independent 
interval approximation (HA). We call this method the Berry- Robnik+IIA 
(BR+IIA) model. In the IIA approximation, the entire statistical behavior 
of the system is described by the nearest neighbor distribution p^ (s) and the 
joint probability to find the particles in positions Xi, ■ • • ,x^ can be written as 
the product of independent product of the N nearest neighbor distributions 



P N (xi, ■ •• ,x N ) = —p {0 \x 2 -x 1 ) ■ ■ •p (0) (x JV -2^_i)p (0) (a;i + L-x j v), ( 33 ) 



lr rhis also can be done for values of k near to 1 or long times by using the numerical 
value of (0) obtained from the simulation and taking into account that (0) = 2q — q 2 
for the Berry-Robnik model. 
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Figure 2: Crossover between the Poisson and Wigner distribution for the 
PCRW, for k = 0.05, for various tirr^s t. The continuous line corresponds 
to the BR+IIA model, the dash line corresponds to the GBR model and the 
square dots correspond to the simulation. 



where the partition function is the normalization constant 

N 

Z N = I dxi . . . dx N J]V 0) (x i+ i - Xi) . (34) 

In the IIA approximation the correlations among intervals [xi, X2], ■ ■ ■ , [xn-i,Xn] 
are neglected, for a review of this method see Ref . [5] . As is natural, the gener- 
alized Berry-Robnik model and the Berry-Robnik+IIA model are equivalents 
in the limit q — > 1 because in this regime those correlations are not strong 
and can be neglected. 

In general, for the PCRW the fits are better for the GBR model than 
the BR+IIA. This means that for this system the information contained in 
p(°'(s) is not enough to describe the entire statistical behavior of the system, 
i.e., the correlation among intervals can not be neglected as it happens in the 
BR+IIA model. These correlations can be neglected only for small times (r < 
ti), when the particles do not coalesce enough to build strong correlations 
between them. The value of T\ as a function of k can be computed from 
equation (13"U1) as follows. Integrating Eq. (13D1) . it is easy to find 

r = 1^* ((S(r)) - (S(0))) + 1 {(S(r)} 2 - (S(O)f) . (35) 

Because of the reaction between particles, (S(r)) grows in time for k > 
making possible to write (S(t)) as (S(0)) plus an increment (AS). In this 
way Eq. ( |35l) takes the form 

r = ] -^- (AS) + - (AS) 2 + - (AS) (S(0)) . (36) 

K 7T 7T 

The correlations between particles can be neglected for small times when 
(AS) 1, then, we can neglect (AS) 2 in Eq. fl36|) . Finally if we consider 
that (AS) must be a fraction e of (S'(O)), we find 

n = e (S(0)) + 1 (Sm 2 \ . (37) 



This equation was derived first in Ref. [9]. From Eq. (l3Tj) . it is clear that 
there is a time where the interaction between particles goes unnoticed even 
in the case of k = 1, where, T\ is the typical time that one particle needs to 
reach one of its nearest neighbors. In the limit k — > 0, we have T\ — > 00 as is 
natural because in this case (S(r)) is a constant. 
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We conclude, that, for the PCRW the spacing distribution and the pair 
correlation functions can be approximated from the uncorrelated superposi- 
tion of a Poisson and a Wigner distribution functions as is proposed in the 
GBR model. 

An alternate way to analyze the crossover of this system, is to study the 
spacing distribution functions at a fixed given time, but for different values of 
k. In order to establish a connection between this picture for the transition 
and the one shown in figure [21 it is necessary to find the correct combination 
of k and r which gives the same statistical behavior for different values of 
these parameters. To find it, in figure [3]-(a) we show the behavior of q as a 
function of r for different values of k. We found that making the change of 
variable 

f = rk 2 /(l - k) 2 , (38) 

all lines shown in figure [3]-(a) collapse in a single one, see figure [3]- (b). Then, 
different combinations of k and r with fixed rk 2 /(l — k) 2 give the same 
statistical behavior, i.e., give the same value for the fit parameter q. This is 
shown in figure HI the data was taken in all cases over 20000 realizations with 
rk 2 /(l — k) 2 « 13.9. For large values of f , f 3> 1, we found that q oc f~ - 5 
and for low values of f, f <C 1, q oc f~ a25 . 

The physical meaning of the change of variable (138]) is clear if we introduce 
the crossover time, T2, between the intermediate regime (which starts when 
t > T\) and the long-time regime (when k renormalizes to 1). This time is 
estimated in Ref. [9], by expanding Eq. (13"TT) in powers of l/y/r, 

T 2 « — p > ( 39 ) 

which is precisely the scaling factor used in figure [3]-(b). The change of 
variable (138]) is f = r/r 2 . 



4 Crossover in the quasi one dimensional gas 
system 

This system was originally studied in [II] . There, the authors studied the 
biased diffusion of two species in a fully periodic 2 x L rectangular lattice 
half filled with two types of particles labeled by their charge, there are L/2 
particles with charge + and L/2 particles with charge — . An infinite external 
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Figure 3: Behavior of the fit parameter q as function of t for different values 
of k. 
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Figure 4: Statistical behavior for the PCRW for different values of k and t 
with rk 2 /(l - kf « 13.9, in all cases q « 0.112. 
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field drives the two species in opposite directions along the x axis (long axis). 
The only interaction between particles is an excluded volume constraint, i.e., 
each lattice site can be occupied by only one particle. The system evolves in 
time according to the following dynamical rules: 

1. L particles are randomly inserted in a 2 x L rectangular lattice, ^ par- 
tides (+) and ^ particles (— ), the remaining sites are empty. Periodic 
boundary conditions are imposed in both directions of the lattice. Let 
the x axis be the long axis of length L. 

2. Two neighbor sites are chosen at random. The contents of the sites are 
exchanged with probability 1 if the neighbor sites are particle-hole, but 
if they are particle-particle the contents are exchanged with probability 
7. The exchanges which result in +/— particles moving in the posi- 
tive/negative x direction are forbidden due to the action of the external 
field. 

3. A time unit, corresponds to 2L attempts of exchange. 

With these dynamical rules, this system evolves with formation of do- 
mains for low values of 7, see figure E3-(b) and E3-(c). In this regime, the 
average length of domains grows in time and while the size of the domains 
remains much smaller than the total size L of the system, the domain size 
distribution exhibits a dynamic scaling. In the long time limit, the sys- 
tem reaches a non-equilibrium steady state (NESS) where there is only one 
macroscopic domain. The length of the macroscopic domain depends on 7, 
for example for 7 = 0.1 it has an approximate size of L/2. Additionally, for 
low values of 7, this macroscopic domain has not a simple charge distribution 
and it almost contains no holes. The macroscopic domain is not in equilib- 
rium because there are particles (travelers) which leak out from one end of 
this domain and travel along the lattice until they reach the other end of this 
domain, see figure [5]- (b). In the case of large values of 7 the system remains 
homogeneous, i.e., disorganized without domain formation, as we can see in 
figure H3-(e). For intermediate values of 7, the macroscopic domain is not well 
formed, it has many holes and it is unstable. In this case there are many 
travelers and small length domains, see figure [5]- (d). 

In order to obtain quantitative results, we measure the length of a do- 
main by using the coarse grained approximation (CG) defined in [TTJ. This 
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I Charge + [J Empty 
(Charge - 



Figure 5: Typical configurations for the gas system in the NESS, for different 
values of 7 from the same initial configuration, (a) Initial configuration, (b) 
7 = 0.1, (c) 7 = 0.3, (d) 7 = 0.5 and (e) 7 = 1 for t = 10000. 
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approximation allows us to map the quasi one-dimensional lattice into a one- 
dimensional lattice. For any configuration on the 2 x L lattice, we construct 
an effective one-dimensional one, with occupation numbers zero or one on a L 
sites line, as follows. At each site i, we assign if there are 5 or less particles 
in the 10 sites around it, including the z-th column of the original lattice. We 
assign 1 otherwise, then, we assign 1 in the i-th site of the one-dimensional 
lattice if there are more particles than holes in the 10 sites around the i-th 
column in the original lattice [E In this simplified description a domain is 
a simple consecutive sequence of ones and its size is just the length of this 
string. 

In figure El we show the behavior of (S(t)) for six values of 7. It is evident 
the NESS behavior in the long time limit when (S(t)) reaches its maximum 
value which depends on 7. Additionally, we found that in the NESS regime 
(S(t)) also depends on the size of the lattice L for all values of 7, in fact 
for 7 = 0.1 it is well know that (S(t)) « L/2. In figured we show (S(t)) 
for 7 = 0.5 and different values of L. The value of (S(t)) increases with 
the value of L and we can expect that it reaches its maximum value in the 
limit L — > 00 . We can conclude that for finite systems (S(t)), in the NESS, 
depends on 7 and L. 

For 7 = 0.1, we found in the scaling regime (S(t)) oc t a6 , this result 
coincides with the one found in [11]. For lower values of 7 it seems that there 
is also a dynamical scaling region whose size decreases when 7 increases. 
Naturally, for 7 = 1 the systems remains homogeneous. Note that (S(t)) is 
very different for 7 = 0.1 and 7 = 0.3 cases, because for 7 = 0.3 there are 
more travelers than in 7 = 0.1, see figure [5) 

In figure [U we show the spacing distribution and the pair correlation 
functions for the gas system for different values of 7. As it can be expected, 
the CG description is not appropriate to measure the length of small domains. 
However, this method allows us to measure of the length of big domains. For 
all values of 7, the nearest neighbor spacing distribution function is well 
fitted by Berry-Robnik model for high values of s. For small values of 7, 
p(°'(s) is well described by the Wigner distribution. In the case of 7 = 1, 
the nearest neighbor distribution is described by the Poisson distribution. 

2 In our description we take only four neighbor lines around the i-th column, two at the 
left and two at the right. This election is arbitrary however the results are not so sensitive 
to it. If we take less neighbor lines we improve the measure of the small domains getting 
worse the measure of big domains and vice versa. Our choice of four neighbor lines allows 
us a reasonable measure of small and big domains. 
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Figure 6: Behavior of (S(t)) for different values of 7, the dashed line is a 
reference line with slope 0.6. We used L = 1000. 

However, the generalized Berry-Robnik model does not describe the next 
spacing distributions nor the pair correlation function with enough precision. 
The differences between the pure coalescing random walk (PCRW for k — 1) 
and the gas system was already studied in Ref. j5] for the case of 7 = 0.1. 
There the authors found that the independent interval approximation (IIA) 
is a better model for the gas system than the CRW model. The results 
that we found for the Berry-Robnik+IIA model are shown in figure [HJ We 
can see that the Berry-Robnik+IIA model is a better approximation for the 
statistical behavior of the gas for all values of 7. This fact suggests that, as 
it happens for 7 = 1, the domains in this system are not strongly correlated 
for all values of 7 and because of that those correlations can be neglected. 

In figure [HI we show the behavior of q as a function of 7 in the interval 
[0.1, 0.35]. We found that at least for 7 > 0.4 the system is in a disorganized 
state, i.e., its statistical behavior is well described by the Poisson distribution 
for large values of s and fits give q ~ 1. In all fits used to find q as function 
of 7 we eliminated the regions where the points are highly dispersed, i.e., 
where the coarse grain method does not measure the length of domains with 
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Figure 7: Average length of domains (S(t)) in the NESS regime for different 
values of L. 

enough precision. 

The numerical results shown in figures [8] and [9] were obtained using a 
2 x 1000 lattice over 20000 realizations, the data was taken at three different 
times t = 2000, t = 3000 and t = 4000. The results shown in figure [6] were 
obtained by using the same lattice over 500 realizations. 

5 Crossover in the spin system 

This system was originally introduced in Ref. [12]. There, the authors con- 
sider a lattice of length L with Lfi spins up ("+") and L(l — ^i) spins down 
("— ") with < /i < 1. Periodic boundary conditions are imposed. The 
spin-flip events are: 

1. + + «->+_+_ A = 4J-£. 

2. --++«->- + -+ A = AJ + E. 

3. + + -+ «-> + - ++ A = —E. 
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Figure 8: Crossover for the gas system. The continuous line corresponds to 
the BR+IIA model, the dash line corresponds to the GBR model and the 
square dots correspond to the simulation. 




Figure 9: Behavior of q in the gas system for different values of 7. 
4. - + «-> +- A = —E. 

The transition probability rate for a process from left to right is for 
A > and 1 for A < 0. The constant J is the nearest neighbor coupling 
between spins, E is the energy associated to an external field which drives 
the up ("+") spins to the right and the down ("— ") spins to the left, and T 
the thermal energy (temperature times Boltzmann constant). 

In Ref. [12] the authors restrict their study to the regime T <C E <C J . In 
this regime the microscopic dynamics of the lattice of spins may be mapped 
into one for an array domain dynamics, which provides a good approxima- 
tion in this regime, for more information see p2, El EJ 121 EE]. With this 
macroscopic description in Ref. [5] the authors show that this system has 
a statistical behavior very similar to the one of the gas system. Addition- 
ally, the system exhibits dynamical scaling behavior and, in the long time 
limit, the system reaches the NESS, where there are two macroscopic do- 
mains which move in opposite directions. However, this domain model does 
not allow us study the crossover of the system between organized and disor- 
ganized states because it is valid only in the regime T <C E <C J- To study 
the crossover regime we must to use the microscopic dynamical rules listed 
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above. In all simulations, we took E = 1 and T = 1. 



(6) 



(d) 



Spin up 
Spin down 



Figure 10: Typical configurations for the spin system in the NESS for differ- 
ent values of J from the same initial configuration, (a) Initial configuration, 
(b) J = 2, (c) J = 1.5, (d) J = 1 and (e) J = 0.1. We use L = 50, // = 0.5 
and t = 100000. 

In figure [TUl-(b). we can see the asymptotic state for the spin system for 
large values of J where there are only two stable macroscopic domains and 
(S(t)) = L/2. For J = 0.1, domain formation is not perceptible, see figure 
[TOl-(e) which it is almost identically to the disorganized initial state figure [101- 
(a). However, as we will see soon, the average length of domains is different in 
both cases. In the NESS, for low values of J there is domain destruction and 
formation, in such way that (S(t)) is constant with a value lower than L/2. 
As we can see in figure [TTl for low values of J, the statistical behavior of the 
spin system in the NESS is well described by the Poisson distribution. That 
means that the system remains in a disorganized state where the average 
length of domains is bigger than the initial one ({S(0)) = 2). We can see 
slight differences near s = 0, this happens because we used in our simulation 
a discrete finite lattice. 

In figure [12], we can see the average length of domains (S(t)) as a function 
of time t. In this case the behavior of (S(t)) is different from the one found 
in the PCRW and in the gas system. For large values of J, the spin system 
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Figure 11: Spacing distribution function in the NESS for J = 1. We took 
data of three different times t/L = 10, 100, 1000 with L = 200 

reaches metastable states where (S(t)) is almost constant. Note that in all 
cases, the metastable state starts when the domain density is 1/4, i.e., when 
the system, in the statistical average sense, is filled with domains of the form 

+ H h H • This is a consequence of the microscopical dynamical rules 

that we use, because they take into account interaction among four neighbors 
spins. For large enough values of J, these domains have a low probability of 
destruction setting the system into a metastable state. In figure [T3l we show 
the spacing distribution function in the metastable states, the level repulsion 
is present. After some time, (S(t)) starts to grow again in time. In the 
PCRW and in the gas systems the metastable states are not present. Once 
the metastable states ends, the domain formation in the system continues 
in such way that the average size of domains increases again. The length in 
time of these metastable states depends on the value of J (with E = T = 1). 
In fact for J = 1 this region is absent but for bigger values of J its size 
increases. Nevertheless, (S(t)) does not seem to depend on the length of 
the system in the scaling regime, it is the same for L = 100, 200, 500, the 
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differences arises near to the NESS naturally H. In figure [T2J for J = 2, 
the system seems to reach another metastable state where (S(t)) ~ 40 for 
N = 200 and (S(t)} « 30 for N = 100. For J = 2, the scaling region has a 
considerable size but this region for lower values of J is smaller. From figure 
[121 it is evident that the value of (S(t)) in the NESS regime depends on the 
size of the system, as it also happens in the gas system. In the inset of figure 
IT2"| we compare (S(t)) for J = and J = 0.1, in the last case there is domain 
formation and it seems that only for J = there is not domain formation. 




Figure 12: Average length for different values of J. Open symbols: L = 100 
and filled symbols: L = 200. In the inset we compare J = and J = 0.1 for 
L = 200 in order to verify the domain formation for small values of J. 

As it happens in the gas system, the statistical behavior of this model 
is better described by the Berry-Robnik+IIA than the generalized Berry- 
Robnik model. In figure [TH we show the statistical behavior of the system 
for different values of J. For J = 0, the system evolves in time remaining dis- 
organized, as we expected. However, we cannot compare directly its spacing 
distributions with our analytical model because of the effects of the discrete 
finite lattice that we use in our simulations. In fact, as we can see in figure 

3 The case L = 500 was verified but is not included in the figure. 
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Figure 13: Statistical behavior of the metastable state in the spin system for 
J = 2. We took data at tree different times: T/L = 10, 20, 30. 

Q3 p(°)(l/2) = 1 instead of p^°\0) = 1 as it is predicted by the continuous 
Poisson distribution. This result can be understood if we remember that, in a 
lattice, the lowest value for which the nearest spacing distribution is defined, 
is s = pit), with p{t) the density of particles. In our case, p(t) = 1/2 for 
J = in each spin species. This technical problem is usually solved taking 
low density values as we did it in the PCRW, where we took p(0) = 1/10. 
However, we cannot do that in the spin system. We have to take small lat- 
tices because the simulation is very intensive and low densities give us a poor 
statistic. This is only a technical problem and we can be sure that in the 
limit of low densities and small values of J, the system is well described by 
the continuous Poisson distribution. To corroborate this, we compare the 
results of our simulation for J = with the discrete version of the Poisson 
distribution (p (0) (s) = (1/2) 25 " 1 , pW(s) = (2s - l)(l/2) 2s ~ 1 , etc.). 

For J = 2 the system sets into an organized state where there is domain 
formation. In this case, p^(s) is well fitted by the Wigner distribution. For 
J = 2, the results of our microscopic simulation coincides with the results 
of the macroscopic array domains simulation in the regime T ^ E ^ J. 
For intermediate values, the system shows a mixed state. If we do not take 
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into account the discrete finite lattice, for all values of J, p( n \s) is well de- 
scribed by the Berry-Robnik+IIA model. This suggest that for this system 
the correlation between domains are not strong as it happen in the gas sys- 
tem. For low densities and values of J near 1, we expect that the generalized 
Berry-Robnik and the Berry-Robnik+IIA models are a good approximation 
for higher spacing distribution functions (n > 1) and for the pair correla- 
tion function because in this regime the correlation between domains can be 
neglected. 




Figure 14: Crossover for the spin system. 
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In our numerical simulations we took /i = 0.5 and 2000 realizations. The 
measures were taken at different times using two lattices: L = 500 and 
L = 200, in order to confirm the scaling property. For J = 2, we took 
t/L = 10000 and t/L = 4000 for L = 500 and L = 200. For J = 1.75, we 
took t/L = 9000 for L = 500 and t/L = 6000, 7500, 9000 for L = 200. For 
J = 1.5, we took t/L = 2000 for L = 500 and t/L = 1500, 2000 for L = 200. 
For J = 0, the system remains disorganized and the time when the data is 
taken is irrelevant. 

6 Conclusion 

In all systems considered in this paper, the original Berry-Robnik model 
gives us a good fit for the nearest neighbor distribution p(°\s), i.e., 
can be approximated by the uncorrelated superposition of one Wigner and 
one Poisson distribution functions. The biggest differences in the fits are 
given when the densities of both sequences have similar values. We calculate 
the next spacing distribution functions by using two methods, the GBR and 
the BR+IIA. We found that the GBR model gives us a good approximation 
for the statistical behavior of the PCRW and the BR+IIA model is a good 
approximation for the gas and spin systems. This fact suggest that the 
correlation between domains in the gas and spin system are not strong and 
can be neglected in a first approximation but in the PCRW the correlations 
between particles cannot be neglected. 

Our analytical models are simple and allow us to have a quantitative mea- 
sure of the degree of order /disorder of the systems through the parameter 
q. For q = 1 the system is completely disordered and p^s) is described by 
the Poisson distribution, that means that the system is homogeneous, in the 
statistical average sense. For q = the system is organized and p^°'(s) is de- 
scribed by the Wigner distribution. In the particular case of the PCRW q = 
implies that there is an effective repulsion between particles and in the case of 
the gas and spin systems means domain formation. Unfortunately, the infor- 
mation about the interaction between particles/domains cannot be extracted 
easily from the statistical behavior in non-equilibrium systems. However, by 
using the BR+IIA model it is possible to calculate in an approximate way 
the interaction potential between domain borders for the gas and spin sys- 
tems, as follows. Comparing (133|) with a Boltzmann factor with an inverse 
temperature (5 — 1, is straightforward to show that, under the BR+IIA ap- 
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proximation, the statistics of the domain borders is approximately equivalent 
to a system with N particles which interact according to the potential 

N N 

V N (xx, ■ ■ ■ , x N ) = q {x i+ i -x i )-} j In [f(xi, x i+1 )] , (40) 

i=l i=l 

where 



f(x t , x l+l ) = g 2 erfc 6^(1 - q)s\ + (2q(l - q) + |(1 - qfs) ',(41) 

and s = — Xi. Thus, the statistics of the domain edges of the gas 
and spin systems in their crossover regimes is approximately equivalent to a 
statistical equilibrium system of particles on a circle interacting through the 
nearest neighbor pair potential PHI) . Note that for q = 0, Eq. (|40p reduces 
to 

N 

V N (xi,--- ,x N ) = ^2^(x i+1 -Xi) 2 -\n(^(x i+1 -Xi)^ , (42) 

i=l 

as we can expect from Ref. [5]. Naturally, for g = 1 the interaction potential 
is a constant and the interaction force between domain edges vanishes. 

In the crossover between order and disorder all systems lose their level re- 
pulsion properties, i.e., the correlations between domains/particles decrease. 
In the PCRW, the correlations between particles arises from the coalescence 
reaction for k > 0, in the gas case, correlations arise from the mutual ob- 
struction of particles for 7 < 1. For the spin system, taking E — 1 and 
T — 1, the crossover depends only on J and the correlations arise for J > 0. 

The gas and spin systems have a similar statistical behavior in the crossover 
between organized and disorganized states but (S(t)) is very different in both 
cases. The metastable regions that we found in the spin system are not 
present in the gas system nor in the PCRW. In fact metastable regions are 
not predicted by the macroscopic dynamical rules used in Ref. [12]. In the 
PCRW there is a time T\ where the interaction goes unnoticed which depends 
on the parameter k. In the gas an spin systems it seems that this time does 
not depend on the parameters 7 and J respectively. 

Finally, for the PCRW, we found that the parameter q = q(r, k), which 
characterize the crossover, can be scaled in a function of a single argument 
f by making the change of variable f = rk 2 /{l — k) 2 . The scaling factor 
is proportional to the crossover time between intermediate regime and long 
time regime, r 2 . 
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